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Abstract — In this paper, we consider the problem of minimizing 
the completion delay for instantly decodable network coding 
(IDNC), in wireless multicast and broadcast scenarios. We are 
interested in this class of network coding due to its numerous 
benefits, such as low decoding delay, low coding and decoding 
complexities and simple receiver requirements. We first extend 
the IDNC graph, which represents all feasible IDNC coding op- 
portunities, to efficiently operate in both multicast and broadcast 
scenarios. We then formulate the minimum completion delay 
problem for IDNC as a stochastic shortest path (SSP) problem. 
Although finding the optimal policy using SSP is intractable, 
we use this formulation to draw the theoretical guidelines for 
the policies that can efficiently reduce the completion delay in 
IDNC. Based on these guidelines, we design a maximum weight 
clique selection algorithm, which can efficiently reduce the IDNC 
completion delay in polynomial time. We also design a quadratic 
time heuristic clique selection algorithm, which can operate in 
real-time applications. Simulation results show that our proposed 
algorithms efficiently reduce the IDNC completion delay com- 
pared to the random and maximum-rate algorithms, and almost 
achieve the global optimal completion delay performance over 
all network codes in broadcast scenarios. 

Index Terms — Network Coding, Instantly Decodable Network 
Codes; Completion Delay; Wireless Multicast; Wireless Broad- 
cast; Index Coding. 

I. Introduction and Motivation 

MULTICAST Broadcast Services (MBS) have become a 
corner stone in the design of all future wireless and 
mobile standards and networks, such as LTE and WiMAX. 
It is now very common to find a group of pedestrians or 
people on a bus watching a soccer match on their smart 
phones while some of them are downloading files, others are 
watching the breaking news and others are looking for the 
closest restaurants or using voice over IP applications,... etc. 
Due to the high demand on these MBS applications and their 
high bandwidth and delay requirements, it is very important 
for MBS protocols to not only efficiently utilize the scarce 
bandwidth resources available to the network but to allow 
progressive packet reception to satisfy delay requirements. In 
other words, while these protocols are minimizing the amount 
of resources (e.g. number of transmissions) consumed by such 
application to increase the bandwidth efficiency, they should 
also be able to guarantee the quality of streaming and delay- 
intolerant applications, in which the received packets should 
be always useful at their reception instant, in order to prevent 
interruption or flickering of the stream. The simultaneous 
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achievement of these goals calls for new approachs to increase 
the efficiency of the packet transmission and recovery pro- 
cesses. One major breakthrough in this area came with the 
development of network coding (NC). 

Recently, NC has shown great abilities to substantially 
improve transmission efficiency, packet recovery, throughput 
and delay over broadcast erasure channels [2]— [16]. Two trends 
of network coding can be distinguished in the literature, 
namely random (or full) network coding (RNC) [17] [18] and 
opportunistic network coding (ONC) [19] [20]. The former 
trend combines all the packets with random non-zero and 
independent coefficients in each transmission. The latter trend 
exploits the receivers' side information in selecting packets 
to be coded in each transmission to achieve a certain target. 
Despite the great interest in RNC in the literature, its ability 
to recover packets without feedback, and its optimality in 
reducing the number of packet transmissions in broadcast 
scenarios [21], it is only feasible for applications with high 
delay tolerance, since it does not support progressive packet 
decoding. It is also inefficient in unicast and multicast sce- 
narios, in which different groups of receivers are interested in 
different subsets of the transmitted packets [14]. 

In this paper, we are interested in a sub-class of oppor- 
tunistic network coding, called the instantly decodable network 
coding (IDNC), in which received packets are allowed to be 
decoded only at their reception instant and cannot be stored for 
future decoding. This IDNC approach is currently attracting 
much attention [8], [9], [22] due to its numerous desirable 
properties. First, IDNC provides instant packet recovery upon 
appropriate packet reception, a property that perfectly matches 
the requirements of the MBS streaming applications, and that 
RNC and general ONC lack. Moreover, the IDNC encoding 
can be implemented using binary XOR, which eliminates 
the complicated operations over large Galois fields and the 
coefficient reporting overhead. This XOR encoding also sim- 
plifies the decoding process at the receivers, and eliminates the 
need for matrix inversion, which represents a computational 
bottleneck in RNC [8]. Finally, no buffers are needed at the 
receivers to store coded packets for future decoding. These 
simple decoding and bufferless properties allow the design 
of simple and cost-efficient receivers, which is an important 
requirement for mobile hand-held devices. 

Despite its attractive properties, IDNC does not guarantee 
service to all receivers in each transmission, which affects its 
ability to both maximize the number of decoding receivers in 
each and every transmission (i.e. capacity) and minimize the 
number of transmissions to deliver a frame or packets (i.e. 
completion delay). The former problem was deeply investi- 
gated in IDNC [8], [9], [22] and ONC [15], [16], whereas 



limited work has addressed the latter problem only for erasure- 
less channels (a.k.a. the index coding problem [2]-[4]). It has 
been shown that this index coding problem is NP-Hard to 
solve and to approximate [23], [24]. Naturally, this complexity 
becomes worse in case of erasure channels, which leaves us 
no choice other than designing efficient heuristics to solve it. 

This fact raises the following question that we address in 
this paper: What is the efficient heuristic policy that can reduce 
the expected completion delay in IDNC, over erasure channels, 
for both multicast and broadcast scenarios? Intuitively, one 
might think that the best heuristic policy is to maximize 
the number of receivers that can decode a new packet in 
each transmission, as studied in [8], [9], [15], [16], [22]. In 
this paper, we show that this intuition is not true and that 
the solution to the completion delay problem is obtained by 
giving priority to targeting the receivers with higher demands 
and worse channels. To reach this result, we first extend the 
IDNC graph, which represents all the feasible IDNC packet 
combinations according to receivers' side information, to suit 
both multicast and broadcast scenarios. We then formulate the 
minimum completion delay problem in IDNC as a stochastic 
shortest path (SSP) problem, which is a special case of the 
Markov decision process (MDP), having absorbing states. 
Although this formulation is impossible to solve, we mainly 
employ it to draw the theoretical properties of the policy that 
can efficiently reduce the completion delay in IDNC. 

Based on these properties, we design a two-stage maximum 
weight clique selection algorithm to reduce the completion 
delay in IDNC in polynomial time for moderate graph sizes. 
For further complexity reduction, we design a quadratic time 
heuristic algorithm, based on greedy maximum weight vertex 
search, which is more suitable for real-time applications. We 
finally compare the performance of our proposed optimal and 
heuristic maximum weight clique search algorithms to RNC, 
the random IDNC algorithm (that selects served receivers ran- 
domly) and the maximum clique IDNC algorithm (maximizing 
the number of decoding receivers). 

The contributions of this papers are summarized as follows: 
• To the best of our knowledge, this paper presents the first 
rigorous study on reducing the IDNC completion delay 
over erasure channels. In this study, we do not limit our- 
selves to maximizing the number of decoding receivers 
as in [8], [9], [15], [16], [22] but rather investigate both 
the order of receiver service and the evolution of coding 
opportunities along the transmission process, which were 
shown to be the key factors affecting the optimization of 
completion delay in IDNC. 
> We design polynomial and quadratic-time heuristics that 

achieve near optimal completion delay performance. 
The rest of the paper is organized as follows. We first 
summarize related works in Section II. In Section III, we 
introduce the system model and parameters. The IDNC graph 
is illustrated in Section IV. We present the problem formula- 
tion in Section V and draw the properties of efficient IDNC 
completion delay reduction in Sections VI, VII and VIII. The 
proposed algorithms are introduced in Section IX and their 
performances are evaluated in Section X. Finally, Section XI 
concludes the paper. 



II. Related Work 

Since its first introduction in [25], network coding has 
been a great attraction to numerous studies as a routing and 
scheduling scheme that attains maximum information flow 
in a network. In [2]-[4], the problem of determining packet 
combinations, to minimize the number of transmissions (i.e. 
completion delay) over erasure-less channels, was studied 
under the name of "index coding". In [23], it has been shown 
that finding the optimal solution of the index coding problem is 
NP-hard and thus different heuristics to solve the index coding 
problem were proposed in [4]. In this paper, we extend the 
study to the case of erasure channels. Our problem differs 
from index coding in that the feedback status of different 
receivers changes probabilistically after each transmission over 
erasure channels. Thus, the coded packets cannot be scheduled 
for the whole transmission process all at once, as in index 
coding, but rather require to be dynamically scheduled after 
each transmission according to the received feedback. 

In [5], the authors proposed an online network coding 
algorithm for the three-receiver case, proved its rate optimality 
and conjectured its asymptotically optimal average delay. 
In [6] and [7], the decoding delay performance of offline 
algorithms was analyzed and the decoding delay of several 
greedy online NC algorithms were compared for i.i.d. erasure 
channels. These proposed algorithms performed un-prioritized 
packet selection for each NC transmission and did not consider 
the channel conditions in their selection procedures. [8], [9] 
proposed a prioritized and channel-aware packet selection 
algorithm that achieves optimal decoding delay for a more 
strict version of IDNC. All these works are clearly different 
from our problem in terms of objective and proposed solutions. 

For a more general ONC scenario than IDNC, in which un- 
decoded packets can be stored for future use, [15], [16] study 
the maximization of service rates (i.e. capacity) of multiple 
unicast sessions over 1-to- -FT broadcast erasure channels. Inner 
and outer capacity bounds were derived and were shown to 
meet in the special cases of symmetric and spatially inde- 
pendent erasure channels. [15], [16] also proposed capacity 
achieving packet evolution algorithms, serving subsets of 
receivers with incremental sizes in order to maximize the 
number of decoding receivers in every transmission, and thus 
their rates. Despite the more general decodability assumption 
considered in these works, our paper differs from them in both 
the optimization objective and proposed algorithms. Unlike 
the aim of [15], [16] to maximize the achievable unicast rates, 
through maximizing the number of decoding receivers in every 
transmission, our paper studies the problem of minimizing 
the total number of transmissions (i.e. completion delay) to 
deliver a frame of multicast or broadcast packets over erasure 
channels, which makes it an extension to the index coding 
problem. Moreover, the proposed packet evolution algorithms 
in [15], [16] do not prioritize receiver service, but rather serve 
receiver subsets, with same size, in an arbitrary sequential 
cyclic or acyclic fashion. Consequently, these algorithms are 
not suitable solvers to the completion delay minimization 
problem, whose solution greatly depends on the prioritization 
of receiver service according to their demand and erasure 
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probabilities, as will be shown in Section VII. Our proposed 
algorithms focus on implementing this prioritization and are 
thus significantly different from those proposed in [15], [16]. 

III. System Model and Parameters 

The system model we consider in this paper consists of a 
wireless sender that is required to deliver a frame (denoted 
by AT) of N source packets to a set (denoted by M) of 
M receivers. Each receiver is interested in receiving either a 
subset or all the packets of M. The former case is referred 
to as "multicast" whereas the latter case is referred to as 
"broadcast". We will refer to the requested and undesired 
packets of any receiver by its "primary" and "secondary" 
packets. The sender initially transmits the N packets of the 
frame uncoded in an initial transmission phase. Each receiver 
listens to all transmitted packets (even the ones that it does not 
want) and feedbacks to the sender a positive acknowledgement 
(ACK) for each received packet. At the end of the initial 
transmission phase, three sets of packets are attributed to each 
receiver i: 

• The Has set (Hi) is defined as the set of primary and 
secondary packets correctly received by receiver i. 

• The Lacks set (d — M \ Hi) is defined as the set of 
primary and secondary packets not received by i. 

• The Wants set (W; C d) is defined as the set of primary 
packets that receiver i has not yet received. 

The sender stores this information in a state feedback matrix 
(SFM) F = [fij] , V i e M,j e Af such that j l3 = if 

j e Hi, fa = l if j e W h and /y = -l if j e d \ W*. 

After the initial transmission phase, a recovery transmission 
phase starts, in which the sender exploits the reception diver- 
sity in the SFM to employ NC. These NC packets must include 
at most one source packet from the Wants or Lacks sets of a 
subset or all of the receivers. The receivers that cannot decode 
a new source packet from this NC packet discard it. For each 
decoded source packet, the receivers send ACK packets that 
are used by the sender to update the SFM and the sets Hi, d 
and Wj, V i. This process is repeated until all receivers obtain 
their requested packets. We define the completion delay of 
a frame as the number of recovery transmissions required to 
deliver all requested packets to their receivers. 

Define g = [g 1: . . . , g M ], ip = [ip 1: . . . , <p M ] and ip = 
[tpi,..., iPm\ as the Has, Lacks and Wants vectors, such 
that Qi, ifi and tpi are the cardinalities of Hi, d and Wi, 
respectively. Let pi and qi = 1 — pi be the packet erasure 
and success probabilities observed by receiver i, respectively. 
We assume that pi and eft V i do not change during the 
frame transmission period. Also, let /i; be the demand ratio 
of receiver i, defined as the ratio of its primary packets in the 
frame to the frame size N. Given this definition of pi, we can 
focus on studying the multicast scenario, since the broadcast 
scenario can be viewed as a special case of the multicast 
scenario, in which /i^ = 1 V i. Finally, define /i = -jjYldLi^i 
as the average of the demand ratios of all receivers. 

IV. IDNC GRAPH 

To form optimized IDNC packets, we should first design 
a representation of all feasible packet combinations that are 



instantly decodable by any subset or all the receivers. An 
initial idea about the representation of packet combinations 
was introduced in the form of a graph, when designing a 
heuristic algorithm to solve the index coding problem [2], [4]. 
This graph, which we will denote by G P (V P , £ P ), is constructed 
by first inducing a vertex Vij in V p for each packet j £ Wi, 

V i £ Ai. Two vertices and Vki in Q p are connected by an 
edge in £ p if one of the following conditions is true: 

> CI: j = 1 => The two vertices are induced by the loss of 
the same packet j by two different receivers i and k. 

• C2: j £ Hk and I £ Hi => The requested packet of each 
vertex is in the Has set of the receiver of the other vertex. 
Consequently, each edge between two vertices in the graph 
represents a coding opportunity, which is defined as an op- 
portunity of generating an instantly decodable packet for the 
two receivers inducing these vertices. Given this graph, we 
can easily define the set of all feasible packet combinations in 
IDNC as the set of packet combinations defined by all maximal 
cliques in Q p (a maximal clique is a clique that is not a subset 
of any larger clique). The sender can generate an IDNC packet 
for a given transmission by XORing all the packets identified 
by the vertices of a selected maximal clique k p in Q p . 

The above formulation of Q p is suitable when optimizing 
packet combinations in a broadcast setting as in [1]. In 
multicast scenarios, we can explore the enhancement of coding 
opportunities at receivers that are not considered for primary 
packet reception, by delivering secondary packets to them. Al- 
though these packets are not requested at these receivers, their 
reception along the steps of the recovery phase, when they 
are not targeted with primary packets, will enlarge their Has 
sets. According to Condition C2, this will increase chances 
of creating more coding opportunities that can serve these 
receivers in the future steps towards completion. However, this 
service of secondary packets should never affect the instant 
decodability of the primary packets at the other receivers. 

To achieve both goals, we propose a new two-layered graph 
Q(V,£). The primary layer consists of graph Q p , described 
above. The secondary layer Q a (y a ,£ a ) is constructed by 
generating a vertex £ V a for each packet j £ C \ Wi, 

V i £ M., and connecting any two vertices satisfying either CI 
or C2. Finally, we connect any two vertices from both layers 
if either CI or C2 holds. In the rest of the paper, we will refer 
to Q p , Q a and Q as the primary, secondary and IDNC graphs, 
respectively. Figure 1 depicts an example of a feedback table 
and its corresponding IDNC graph. It is easy to show that the 
overall complexity of graph construction is 0(M 2 N). 

According to the design of Q, we can easily infer that each 
receiver can have at most one primary or secondary vertex 
in each of its maximal cliques. Consequently, the selection 
of a maximal clique for a given transmission is equivalent 
to the selection of a set of targeted receivers. We can thus 
envision the completion delay minimization problem, given the 
IDNC graph, as a problem of finding the optimal scheduling of 
targeted receiver subsets under the IDNC constraint, in order to 
complete the recovery phase as fast as possible. After selecting 
the maximal clique determining the targeted receivers for any 
transmission, the source packets identified by its vertices will 
be XORed to generate the coded packet. 



3 



Packet 1 Packet 2 Packet 3 Packet 4 



Receiver 1 


1 


-1 


1 





Receiver2 





1 


-1 


1 


Receiver3 


1 








-1 












t 


11 J 


f \\ 

' 13 


i f 12 




( 31 1^. 








^{ 34 




[ 22 ) 


(24 \ 

\ J * 
x. y / 


if 23 J 





Primary 
Graph 



Secondary 
Graph 



Fig. 1: Example of a feedback matrix and its corresponding 
IDNC multicast graph. The shaded and white boxes and ver- 
tices represent the wanted and unwanted packets, respectively. 



In the rest of the paper, we say that a receiver is targeted 
by an IDNC packet if the selected maximal clique includes a 
vertex induced by this receiver. We also define T p {n), T a (n) 
and T(k) as the set of all primary, secondary and overall 
targeted receivers of a given maximal clique k. 

V. Problem Formulation using SSP 

A. The SSP Problem 

The stochastic shortest path (SSP) problem is a special 
case of the infinite horizon MDP, which can model decision 
based stochastic dynamic systems with terminating situations. 
In SSP, all the possible situations the system may encounter are 
modeled as states in a state space <S. In each state s G S, the 
system must select an action a from an action space A(s) C A 
that will charge it an immediate cost c(s, a) (where A denotes 
the action space of SSP). The terminating situations of the 
system can be thus represented as zero-cost absorbing states. 
Once an action a is taken at state s, the system can move to a 
state s' with probability P a {s, s'), which only depends on the 
current state and the taken action. An SSP policy tt = [tt(s)] 
is a mapping from S^i that specifies a given action to each 
of the states. The optimal policy tt* of an SSP is the one that 
minimizes the cumulative mean cost until an absorbing state 
is reached. 

The algorithms that solve SSPs define a value function 
K-(s) as the expected cumulative cost until absorption, when 
the system starts at state s and follows policy tt. It can be 
recursively expressed V s 6 S as: 

V T (S) = C(«,7T(«)) + ^(s)MK(s') , (1) 

where S(s, a) is the set of successor states to s when action a 
is taken (i.e. S(s, a) = {s'\P a (s, s') > 0}). Consequently, the 
optimal policy at state s can be defined V s 6 S as: 




arg mm 



s 1 £S(s,a) 



(2) 



Fig. 2: State representation, action space of the example in 
Figure 1 and its possible transitions for action 07. 



B. Problem Formulation 

The problem of minimizing the expected completion delay 
in IDNC can be formulated as an SSP problem as follows: 

1 ) State Space S: States are defined by all possibilities of 
SFM F(s) that may occur during the recovery transmission 
phase. For state s, the matrix represents the content of Has, 
Lacks and Wants sets in s (i.e. Hi(s), Ci(s) and W;(s) 
Vie M) as defined by (??). According to its definition, the 
state space has a size of \S\ = O (2 AIN ). We can characterize 
each state s by its Has, Lacks and Wants vectors q(s), 
<p(s) and ip(s). Note that several states can have the same 
cardinality vectors. The Wants vector of any absorbing state 
is 1 x AI, which we denote by i/j . 

2) Action Spaces A(s): For each state s, the action space 
A(s) consists of the set of all possible maximal cliques in 
the IDNC graph Q(s), constructed from the SFM F(s). 

3) State-Action Transition Probabilities: To define the 
state-action transition probability P K ( s )(s,s') for an action 
k(s) £ A(s), we first introduce the following two sets: 

X = {1 e T(k(s)) I ip^s) > <pi(s')} (3) 

y = {i e T(k(s)) \ <pi(s) = <pi(*')} ■ W 

The first set includes the targeted receivers whose Lacks 
sets have decreased from state s to state s', and thus have 
successfully received the IDNC packet generated from k(s). 
The second set includes the targeted receivers that have lost 
the IDNC packet generated from k(s) and thus their Lacks 
sets did not change. Based on the definitions of these sets, 
-P«(s)( s ) s ') can be expressed as follows: 

P K ( s )(s,s') = IJ qi J} Pi . (5) 

Figure 2 depicts the state representation and the action space 
for the example in Figure 1. It also depicts the possible 
transitions given that action 0,7 is performed. 

4) State-Action Costs: The expected completion delay is 
defined in SSP terms as the expected number of transitions in 
the process before arriving to an absorbing state. Since any 
transition (due to any action) takes one packet transmission, 
the cost payed by the process is one time-slot. Consequently, 
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the costs of all actions in all states should be set to 1 . In other 
words, c(s, k(s)) = 1 V k(s) e A(s), s e S. 

C. SSP Solution Complexity 

The optimal policy of an SSP problem can be computed 
using the policy iteration and value iteration algorithms. The 
complexities of these algorithms are (|5| 3 + |<S| 2 |„4|) and 
(|<S| 2 |„4|). According to the dimensions of S and A(s) 
described in Section V-B, we conclude that computing the 
optimal policy is very difficult in real-time for typical values 
of M and N. Even the simulation based technique proposed 
in [11] will not be able to compute the optimal policy in real- 
time since its complexity still scales with \S\. 

VI. SSP Properties 

Despite the complexity of solving the SSP problem formu- 
lated in Section V-B, we can study its properties and structure 
to draw the characteristics of policies that can efficiently 
minimize the expected completion delay. From Section V-B, 
it is easy to infer that the SSP formulation has the following 
properties: 

Property 1 (Uniform Cost). 

The costs of all actions in all states are all the same except 
for the absorbing state. 

Property 2 (Non-singleton acyclicity). 

No state can be revisited once the process moves to a next 
state, and thus the SSP formulation is acyclic. Indeed, if some 
packets are received by some receivers when an action is taken 
at a given state, there is no means of going back with these 
receivers not having these packets. However, a state can revisit 
itself (singleton cycles) if none of the targeted receivers by the 
taken action receives the IDNC packet. 

Property 3 (Non-increasing successor value functions). 

Since there are no cycles of size more than one, the successor 
states of a state s are all closer to the absorbing states than 
s. Consequently, the expected cost to absorption starting from 
a given state is always greater than or equal to the expected 
costs to absorption starting from all its successor states. 

These three properties can be employed to draw the proper- 
ties of the optimal policy it* minimizing the mean completion 
delay at any given state s as follows. From the uniform cost 
in Property 1, we have: 

tt* (s) = arg min J 1 + V P K {s, s') V** (s') 



irg mm < V P K (s, s') T4» (s') 
y s' eS(s,K(s)) 

k(s)S-4(s) 



(6) 



where E K ( S .) is the expectation operator over the different 
transition probabilities when action k(s) is taken. Thus, the 
optimal action at state s is the action minimizing the expec- 
tation of the optimal value functions of the successor states. 



From Properties 2 and 3, we know that all successors of state 
s are closer to the absorbing state (thus having smaller mean 
completion delays) except for itself. Consequently, the optimal 
action at state s is the one that has high probability in moving 
to states with the minimum expected residual completion delay 
(i.e. minimum mean time to absorption), given the optimal 
policy. 

Now the problem is that there is no close form expression 
for the optimal value functions V n *(s') in IDNC and thus 
there is no means of accurately computing it to determine the 
optimal policy without solving the SSP. However, based on the 
previous properties and facts, we can easily infer that the value 
of V-n* (s 1 ) for any s', that is successor to state s, depends on 
two main factors: 

• The closeness of the state's Wants vector il>(s') to that 
of the absorbing states ipa. 

• The number and sizes of the primary maximal cliques 
available as actions in state s'. 

Indeed, the smaller the distance between if)(s') and i/jo, the 
smaller the value function of state s'. However, this condition 
is not enough as we should also check the availability of 
efficient actions at this state that can bring the system faster to 
an absorbing state. In general, the successor states of s, whose 
primary graphs include more numerous and larger maximal 
cliques, have more chances of reaching the absorbing state 
faster than the others. Since all states s' are successors of a 
same state s, their graphs are different variants of Q(s) depend- 
ing on the vertices that have been served. Consequently, the 
action at state s, which can maximize the coding opportunities 
in the IDNC primary graph at state s', will result in larger and 
more numerous primary maximal cliques, which will help in 
reaching an absorbing state faster in future steps. Based on 
these observation, we state that the policy that can efficiently 
reduce the expected completion delay in IDNC should aim, at 
any visited state, to both: 

« Bring the system Wants vector the closest to the absorb- 
ing states vector i/jo. 

• Maximize the coding opportunities in the successor 
state's primary graph. 

If we can find a policy that can simultaneously achieve these 
two goals, we will employ it to design an algorithm to 
efficiently reduce the expected completion delay for IDNC. 
To investigate the existence of such policy, we will study 
two important features of the problem, namely its geometric 
structure, and the evolution of coding opportunities in the 
IDNC graph. This will be the target of the next two sections. 

VII. Geometric Structure 

In this section, we will explore the actions, which have high 
chances of moving the system Wants vector closest to that of 
the absorbing states. Given the representation of the SSP states 
by their Wants vector, we can define a geometric structure as 
follows. Define an Af-dimensional space, and locate to each 
point ^ = [^i, . . . , m this space all the states having 
Wants vectors equal to the coordinates of this point. Although 
many states can share the same Wants vector and thus can 
be located at the same point, these states differ from one 
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Fig. 3: Geometric structure of the example in Figure 2 after 
removing the fourth column 



another by their IDNC graphs. All absorbing states will be 
located at the origin ^ of this space. Note that this geometric 
representation has the same non-singleton acyclicity property 
as the SSP (i.e. a point cannot be revisited after it is left). 

Since at most one packet can be decoded by each receiver 
from any IDNC transmission, the system can at most move 
from the point = i/)(s) to a point ty' = ip(s') which is a 
vertex in the hypercube T(s) defined as: 



r(s) = I ^i(s) - Vz(s') G {0, 1} V * e M 



(7) 



In other words, T(s) is the hypercube of side length 1, in which 
ipi(s) i s the corner having the largest coordinates. In this case, 
the optimal action at any state is the one that can transitions 
the system to the opposite diagonal point in the M-dimensional 
hypercube, for which ipi(s) — ipi(s') = 1 V i € M. This 
action means that all M receivers are targeted with primary 
packets. If such actions exist and are applied in all visited 
states, we will reach completion faster. However, these actions 
will most probably not exist in most states due to the instantly 
decodability constraint. Consequently, We need a method to 
estimate the closeness of other points to the absorbing state to 
evaluate our IDNC scheme. 

Figure 3 depicts the geometric structure of the example 
in Figure 2 after removing the fourth column (i.e. removing 
the fourth packet and the actions it appears in). Conse- 
quently the system is at point identified by the Wants vector 
ip(s) = [2,1,1]. In this example, there are only five actions 
ai, <22, d3, as, a§ according to their notation in Figure 2. As- 
suming that the system is erasure-free, action (13 will lead 
the system to point [1,1,1] whereas action 05 will lead it 
to point [2,0,0]. Although 05 targets more receivers than 03 
(which makes it a capacity achieving action as of [15], [16]), 
we can clearly see that the latter action gives the chance 
to the system to reach absorption with one more erasure- 
free transmission, if there exists an IDNC packet targeting 
all three receivers at the state located at point [1,1,1]. This 
closeness to absorption in terms of Wants vectors is shown 
through the smaller geometric distance from point [1,1,1] 
to compared to point [2,0,0]. We can infer from this 
example that minimizing the maximum entry of the Wants 
vector (i.e. m&Xi{if}i(s')}) brings the system closest to the 
absorbing point. The intuition behind this finding is that the 
receivers having the largest Wants sets will impose their Wants 



set cardinalities as lower bounds on the completion delay. 
Consequently, serving these receivers first gives hope to reduce 
this lower bound at each step whereas ignoring them will not 
change the lower bound. 

However, minimizing the maximum of the Wants vector 
entries is not enough to describe the actions with closest 
successor states to absorption. For example, actions a\ and 
03 in Figure 3 have the same value for m&Xi{ipi(s')}, but a\ 
brings the system closer to the absorbing point compared to 03 
in terms of Wants vector, since it serves an additional receiver 
with smaller Wants set. This is also reflected on the geometric 
distance from the two destination points to ^q- 

From the above example, we can conclude that, in order 
to bring the system closest to the absorbing point in terms of 
Wants vector, the sender should give more weight to serving 
the receivers with largest coordinate entries, while maximizing 
the number of served receivers with the smaller coordinate 
entries. This weighting can be done through norm expressions. 
For example, the L2 norm (Euclidian distance) represented the 
state closeness to absorption, in terms of the Wants vectors, 
in the previous examples. The larger the employed norm, the 
more biased the weighting in giving service to the receivers 
with largest coordinate entries. 

For erasure channels, the effect of packet erasures should be 
reflected on the geometric structure of the problem. Let i and 
k be two receivers having the same Wants set size but pi > 
Pk- Consequently, i will require on average more targeting 
attempts compared to k in order to deplete its Wants set. Since 
we assume that erasure probabilities do not change during the 
transmission of a frame, targeting k and ignoring i is expected 
to result in a higher overall completion delay, especially when 
ipi(s) is among the largest value in ij)(s). According to these 
facts and the intuition explained above for the erasure-free 
case, i should be given a higher priority of service than k. 

To implement the above prioritization, we define a channel 
weighted Wants vector ip(s) = tpx(s), . . . ,i/jm(s) , where 

= g ■ Based on this new vector definition, we can 
redefine our space such that its points "J are identified by the 
coordinates of the vectors t/>(s) instead of tp(s) V s G S. In 
this case, the actions move the system within hyper-rectangles 
r'(s) with sides equal to q~[ in the i-th dimension. In other 
words: 



r'( s ) = |*' 



(8) 



The sender should then take the action that can reach successor 
states with minimum L n norm over this geometry. 

From the above observations, we can draw a conclusion that 
the policies, which efficiently reduce the IDNC completion 
delay, should always aim at each visited state s to reach a 
state s' that is located at the point with minimum distance to 
the absorbing point Vfo an d thus the minimum ipi(s') 

n 

Consequently, the receivers with larger values of ijji will 
have higher priority to be selected for transmission at state 
s. Now, if we can show that this norm based selection of 
the receivers also maximizes the coding opportunities in the 
successor states, then this norm based selection policy is our 
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searched policy as explained in Section VI. To investigate this 
point, we will study the evolution of the coding opportunities 
in the IDNC graph with respect to the selection of receivers 
in each transmission in the next section. 



VIII. Evolution of Coding Opportunities 

As stated in Section VI, one major factor that identifies 
the efficiency of an action in reducing the completion delay 
is its ability to maximize the coding opportunities in the 
primary IDNC graph. We know from the IDNC graph structure 
that the coding opportunities are represented by its edges. 
Consequently, the overall number of coding opportunities in 
the graph are represented by its edge set size. In [26], we 
studied the receiver selection strategies, maximizing the coding 
opportunities and density in the IDNC graph for wireless 
broadcast. This study is done based on ignoring the packet 
content of the Has, Lacks and Wants sets and computing the 
expected edge set size given their cardinalities. In this paper, 
we will extend this study to the case of wireless multicast 
by first introducing the following theorem. In all upcoming 
analysis in this section, we assume the system is at an arbitrary 
state s and all the variables represent their values at this state, 
thus dropping the (s) notation for simplicity. 

Theorem 1. For a given maximal clique k, chosen for 
transmission at time t in a multicast scenario, the expected 
edge set cardinality of the IDNC primary graph at time t + 1 
is expressed as: 



E 



= E 
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(0 



at time t. 



is the expected degree of a vertex of receiver i 



Proof: The proof can be found in Appendix A. ■ 
From the above theorem, we can draw the following obser- 
vations about the maximization of the edge set size. 



A. Causal Evolution 

The first term in the right-hand side of (9) is the expected 
edge set cardinality at time t. This means that the edge set 
size at one time depends on its previous size instances and 
thus the evolution process is causal. Consequently, if the edge 
set cardinality at t + 1 was maximized, the future evolutions 
in visited successor states at times t + 2 and t + 3, ... will also 
benefit from the maximization that occurred at time t, which 
results in the maximization of their edge set size, if the same 
policy is employed in each step. 

B. Vertex Disappearance 

The second term in (9) represents an expected reduction 
in the edge set size due to the possible disappearance of the 
primary targeted vertices. This disappearance results in the 
removal of their adjacent edges at time t, which is reflected 
in the E [A*] term. It also results in the loss of the potential 
improvement in these degrees if they were kept in the graph, 
which is reflected in the 7, term. This loss is a natural outcome 
of the recovery transmission process and is unavoidable. We 
cannot try to reduce this term by reducing the size of the 
primary targeted receiver set as this will tend to increase the 
expected completion delay. However, we can still reduce the 
effect of this loss component by serving the vertices with 
smaller degrees. The following theorem compares the expected 
vertex degrees of two receivers given the sizes of their Has 
and Wants sets. 

Theorem 2. If i> l > 4>h an d Qi < Qh, then E [A h ] > E [A 4 ]. 

Proof: The proof can be found in Appendix B ■ 

Now, if q 4 < q h and ^ > iph, \<H^ [\] < §<Z/*E [A ft ]. 
Consequently, serving receivers with largest Wants sets and 
erasure probabilities, and smallest Has sets, results in a smaller 
loss in the resulting edge set size. 

C. Degrees of Remaining Vertices 

The third and fourth terms in (9) represent the change in 
the degrees of the remaining vertices in the primary graph, 
which are quantified by ipicti and ipiPi f° r tne targeted and 
non-targeted receivers, respectively. The following theorem 
describes the relation between these two terms. 

Theorem 3. The increase in the degrees of the remaining 
vertices of any receiver is larger when it is targeted than when 
it is not. In other words, on > f3i V i € M.. 

Proof: The proof can be found in Appendix C. ■ 
Now, moving a receiver i from the non-targeted set to 
the targeted set results in an extra number of edges of 
4r {o-i This term is larger when is larger, and thus 

moving a receiver with a larger Wants set to the targeted 
receiver set adds more edges to the primary graph than moving 
a receiver with a smaller Wants set. Consequently, a larger 
increase in the expected edge set size is obtained when 
targeting the maximum number of receivers having larger 
Wants sets. 

Another important insight about the values of a t and /?, can 
be inferred from the analysis of their components $^ and A^. 
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Since the terms £fceT P (K) $ifcfe) and Y,keT P {n) $ik(0) are 

fc^i k^i 

sub tractive terms from 014 and respectively, then selecting 
the receivers with smaller values of &ik(qi) and $jfc(0) to be 
primary targeted receivers increases the values of a% and Pi, 
respectively. Now, if q k < q h , ijj k > tp h and g k < g h , we 
have: 

qk (Qk - ipk. + 1) < qh (Qh - ^h + 1) (15) 
=> < and $ lfc (0) < $ ifc (0) (16) 

Consequently, the receivers having larger Wants sets and 
erasure probabilities, and smaller Has sets, have smaller values 
of ^ik(q-i) and $jfe(0). In case of equal demand ratios at all 
receivers (including the broadcast scenarios), the three above 
conditions are equivalent on average. In other words, the 
receivers having smaller reception probabilities will on average 
have larger Wants sets and smaller Has sets. Consequently, 
maximizing the number of such receivers in the set of primary 
targeted receiver maximizes the values of both cti and /?,;. 
Since the terms EfceT,(«) A lk ( qi ) and EfceT^) A tk (0) 

k^i k^i 

are additive terms to c*j and j3i, respectively, then selecting 
the receivers with larger values of Aik(qi) and Ai k (0) to be 
secondary targeted receivers increases the values of a* and Pi, 
respectively. The values of A^qA and Ajfc(O) are larger for 
receivers having larger value of q k ip k and thus having larger 
Wants sets but lower erasure probabilities. 

D. Overall Maximization Strategy 

From the above theorems and discussion, we can infer that, 
after a given transmission, the value of the edge set size in 
the primary graph is maximized by targeting the receivers, 
having the largest Wants sets and erasure probabilities (i.e. 
largest ipi values), with primary packets, then targeting the 
remaining receivers, having the largest Want sets and lower 
erasure probabilities, with secondary packets. We will refer to 
this strategy as the worst receiver layered targeting (WoRLT) 
strategy. 

In the primary graph, the WoRLT strategy is equivalent to 
a norm minimization of the channel weighted Wants vector 
i/\ Indeed, such minimization will result in targeting the 
maximum number of receivers having the largest Wants sets 
and erasure probabilities. According to the discussion in 
Section VII, this policy perfectly matches the policy bringing 
the system the closest to the absorbing point ^q. 

For the secondary graph, the WoRLT strategy is still equiv- 
alent to a norm minimization of the channel weighted Wants 
vector i/> within the IDNC secondary subgraph, adjacent to 
all selected vertices in the primary graph. Since the receivers 
with the largest Wants sets and erasure probabilities will be 
targeted with vertices from the primary graph, and since each 
receiver can have at most one vertex per clique, applying 
the norm minimization in the secondary subgraph, adjacent 
to the primary selected vertices, will result in targeting the 
remaining receivers with largest Wants sets and lower erasure 
probabilities. This step does not conflict with the policy 
bringing the process closest to absorption but rather fosters it 
in future steps. Indeed, serving these receivers with secondary 



packets will increase the coding opportunities of their remain- 
ing primary (and secondary) vertices, which allows them to get 
served faster in future transmissions, thus bringing the process 
closest to absorption. 

Given the above facts, we conclude this section by stating 
that the WoRLT strategy is the policy that can efficiently 
reduce the IDNC completion delay, as claimed in Section VI. 
We will thus design our proposed algorithms according to this 
strategy in the following section. 

IX. Proposed Algorithms 

A. Maximum Weight Clique Selection Algorithm 

According to the findings of the previous sections, we 
propose a two-step maximal clique selection algorithm that 
should be executed at any visited state s. In the first step, 
the algorithm selects the maximal clique k*(s) in the primary 
graph that targets receivers with larger channel weighted Want 
set sizes, thus minimizing ||"0i(s)||ri for the expected successor 
state and maximizing the number of edges in its graph. To 
further maximize the number of edges, the same process 
should be done for the secondary subgraph connected to k* 
to find k*(s). Each of these two steps can be done using a 
maximum weight clique selection algorithm as follows. 

For each vertex Vij in the multicast IDNC graph, we assign a 

weight (^pi(s)j , where n is the order of the selected norm. 
After finding this clique, the secondary subgraph connected 
to k* is extracted and the maximum weight clique selection 
algorithm is run on it to obtain «*. When both cliques are 
found, the sender sends an IDNC packet that is generated by 
XORing all the source packets identified by the vertices in both 
cliques. After receiving the feedback from the receiver, the 
sender determines the reached successor state and the whole 
procedure is re-executed. This loop is run until all vertices in 
the primary graph are depleted. 

It is well known that the maximum weight clique selection 
problem is NP-hard [27], and is hard to approximate [28]. 
On the other hand, there exist several algorithms that solve 
this problem in polynomial time for moderate size graphs 
( [29] and references therein). However, the complexity of 
these algorithms may still be prohibitive for the applications 
of interest in this paper [29]. Consequently, we will design a 
simple heuristic in the next section to solve the problem with 
much lower complexity. 

B. Maximum Weight Vertex Search Algorithm 

In this section, we design a simple algorithm that performs 
clique selection, using a maximum weight vertex search. For 
this search to be efficient, the vertices' weights must not only 
reflect the (^ipi(s)j values of their inducing receivers, but 

also their adjacency to vertices having high (^i(s)j . 

To design the vertices' weights, we first define ay.jwfs) as 
the adjacency indicator of vertices and v k i in G(s) such 
that: 



CLij,kl(s) 



1 is connected to Vu in Q(s) 







u %3 

otherwise . 



(17) 



We then define the weighted degree A.y(s) of vertex Vij as: 

**M= E *M*) (Ms))" ■ ( 18 > 

Thus, a large weighted vertex degree reflects its adjacency to 
a large number of vertices belonging to receivers with large 
values of (^tpi (s)J . We finally define the vertex weight Wij (s) 

Wij (s) = Ag(s) . (19) 

Consequently, a vertex u.y has a large weight when it both 
belongs to a receiver with large Ctpi (s) J value and is adjacent 

to a large number of vertices with large (^ijjk(s)j values. 

Based on these definitions, we can introduce our proposed 
packet selection algorithm as follows. The algorithm operates 
only for visited states. In each visited state s, the algorithm 
first computes a primary maximal clique k*(s) in Q p {s). At 
first, k*(s) and «*(s) are empty sets. The algorithm starts 
by selecting the maximum weight vertex in Q p (s) to be the 
source vertex in n*(s). For each of the following iterations, 
the algorithm first recomputes the new vertex weights within 
the primary subgraph connected to all previously selected 
vertices in k*(s), then adds the new maximum weight vertex 
to it. When there is no further primary vertices adjacent 
to all vertices in K*(s), the same process is repeated with 
the secondary subgraph adjacent to k* until no vertices are 
remaining in the global graph. The final maximal clique k* (s) 
is thus the union of k*(s) and ft*(s). Once this clique is 
computed, the sender forms and sends an IDNC packet by 
XORing the source packets identified by the vertices in k*(s). 
According to the received feedback, a new state is visited 
and the process is re-executed until the absorbing state is 
reached. Since a maximal clique can have at most M vertices 
(we can target each receiver at most once per transmission), 
and since each iteration in the algorithm requires weight 
computations for the O(MN) graph vertices, the complexity 
of the algorithm is 0(M 2 N). 

X. Simulation Results 

In this section, we present simulation results comparing 
the performance of our proposed algorithms to the following 
algorithms in both multicast and broadcast scenarios: 

• Random clique search algorithm (RND), employed in 
[13], which picks a random clique from the graph for 
each transmission . 

• Maximum clique selection algorithm (MC), which selects 
a primary maximum clique K 1 p nax (s) in Q p then selects the 
secondary maximum clique from the secondary subgraph 
connected to K™ aa: (s) (denoted by MC in figures). 

• Prefect RNC, in which we assume full independence 
between all transmitted coding coefficient vectors. Thus, 
this scheme represents the global optimal completion 
delay in the broadcast scenario. 

For our proposed algorithm, we consider the L\, L2, L3, L$ 
and Lio norms to test the effect of the selection bias on the 
algorithms' performance. In our simulations, we assume that 
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Fig. 4: Performance comparison of different schemes vs \i 
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Fig. 5: Performance comparison of different schemes vs M 

different receivers have different packet erasure probabilities 
and different demand ratios, that change from frame to frame 
while keeping the average erasure probability (p) and average 
demand ratio (/i) constant. 

Figures 4 5, 6 and 7 depict the average completion delay 
performance of the maximum weight clique selection algo- 
rithm with different norms and compare it to the different 
algorithms against the average demand ratio // (for M = 60, 
N = 30, p = 0.15), the number of receivers M (for \i = 0.5 
and 1, N = 30, p = 0.15), the number of packets N (for 
/i = 0.5 and 1, AI = 60, p = 0.15), and the average erasure 
probability p (for fx = 0.5 and 1, M = 60, N = 30), 
respectively. 

From all these figures, we can draw the following observa- 
tions: 

• Our proposed maximum weight clique selection algo- 
rithm with all norms considerably outperforms the RND 
and MC selection algorithms in terms of average com- 
pletion delay for all comparison parameters (p, M, N 
and p), Note that the MC is the IDNC version of 
the algorithms proposed in [15], [16], which shows the 
superior performance of our algorithms in the context of 
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Fig. 6: Performance comparison of different schemes vs N 



Fig. 8: Performance comparison of optimal and heuristic clique 
search algorithms vs /i 
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Fig. 7: Performance comparison of different schemes vs p 



completion delay. 

• The Li norm algorithm employed in [1] degrades in 
performance, compared to the other norms, since it rep- 
resents a very loose upper-bound of V v *{s). 

• For norms higher than 1, the algorithm tends to converge 
to the same performance with the smallest completion 
delays achieved by the L3 and L5 algorithms. For greater 
norms like L\q, the performance slightly degrades for all 
comparison parameters. 

« For the broadcast case (/1 = 1), results show that our 
proposed algorithm almost achieves the optimal perfor- 
mance of random network coding for all comparison 
parameters, with a maximum degradation of less than 
5% only occurring at very high number of receivers. 
This near-optimal performance is achieved while fully 
preserving the benefits of IDNC compared to perfect 
PvNC. 

Figures 8, 9, 10 and 1 1 compare the average completion de- 
lay of our proposed optimal maximum weight clique selection 
(denoted by "opt") to that of our proposed maximum weight 
vertex search (denoted by "srh") algorithm for i 3 , L 5 and 
Lio, as well as the maximum clique algorithm, using the same 
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Fig. 9: Performance comparison of optimal and heuristic clique 
search algorithms vs M 



simulation parameters in Figures 4 5,6 and 7, respectively. For 
the MC approach, the heuristic algorithm is the same as the 
one described in Section IX-B, in which the \M>ij value of 
a vertex is replaced by its absolute primary degree. 

For all these figures, we can see that the heuristic algorithms 
perform very closely to the optimal clique selection algorithms 
for all norms and all comparison parameters, with a maximum 
degradation of less than 5% only occurring for large numbers 
of receivers. We can also observe a considerable improvement 
of our proposed heuristic algorithm with all norms compared 
to both the optimal and heuristic MC algorithms. 

XI. Conclusion 

In this paper, we studied the problem of minimizing the 
completion delay in wireless multicast and broadcast scenarios 
for IDNC. We formulated the minimum completion delay 
problem for IDNC as a stochastic shortest path (SSP) prob- 
lem and showed that it is intractable. Nonetheless, we were 
able to draw the theoretical properties of the policies, which 
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Theorem 4. For given q, ip and xp vectors, the expected edge 
set cardinality of the primary graph is equal to: 



E[|£, 
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Theorem 5. For a given maximal clique k, chosen for 
transmission at time t, the expected primary degree of a vertex 
of receiver i S T(k) at time t + 1 is expressed as: 
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Fig. 11: Performance comparison of optimal and heuristic 
clique search algorithms vs p 



can efficiently reduce the IDNC completion delay, using the 
properties of the formulated SSP and the nature of the IDNC 
graph evolution. Based on these properties, we designed an 
efficient IDNC algorithm using a maximum weight clique 
selection algorithm, which can be solved in polynomial time. 
For further complexity reduction, we also designed a quadratic 
time heuristic algorithm, which can operate in real-time appli- 
cations. Simulation results showed that our proposed heuristic 
can both efficiently track the optimal performance and outper- 
form the random and maximum clique algorithms. For the 
broadcast case, simulations also showed that our proposed 
algorithms almost achieve the globally optimal performance 
of random network coding, while preserving all the benefits 
and simplicity of IDNC. 

Appendix A 
Proof of Theorem 1 

To prove this theorem, we first introduce the following two 
theorems, proved in Appendices D and E, respectively. 



When the maximal clique k is chosen for transmission at 
time t, each member k of the targeted receiver set T(k) may 
(may not) receive the coded packet with probability qk {pk)- 
Let Xk be the random variable representing the reception of 
receiver k G T(k) at time t and X as the random vector of 
all such random variables. From Theorem 4, we can derive 
the expression of the expected edge set size at time t + 1, 
conditioned on the random vector X, as follows: 
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Now, taking the expectation operator over the random vector 
X, and using the expressions of Theorem 5, we get the 
expression of the expected edge set size at time t+1 as follows: 
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Grouping similar terms and expanding E 
get: 
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Appendix B 
Proof of Theorem 2 

To prove this theorem, we first introduce the following 
theorem, proved in Appendix F. 

Theorem 6. For given q, (p and if) vectors, the expected 
primary degree of any of the vertices induced by receiver i 
(denoted by AJ is equal to: 
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Expanding the expression in (25), we get: 
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Appendix C 
Proof of Theorem 3 



_\p qii>kQk \p gfcV'fc (a + gi) 
k =£ % k^i 

Ik A . (gfc - ^fc + 1) (ft + 9») 



keT P (K) 

k^i 



N - 1 



(27) 



Re-arranging the above expression and using the definition of 



(Since for ip^ > 0, Qf. > Qk — tpk + 1 an d V'fc > 9fc> tne l ast 
'term in (28) is non-negative and the theorem follows. 

Appendix D 
Proof of Theorem 4 

It is well known from graph theory that the edge set size 
of any graph is equal to half the sum of its vertex degrees. 
Consequently, we will find an expression for the sum of the 
vertex primary degrees as follows. From Appendix F, we know 
that the primary degree of a vertex Vij is expressed as in (37). 
Consequently, the sum EA^ of all the primary degrees of the 
vertices induced by receiver i can be expressed as: 

M 

SA > = E E 



M 

E 

k=l 

k^i 

M 

E 

fe=i 

k^i 



E ^ ew <= 



I jeWk +I je K k (\Wk\-\mn£i\) 

- Ijen k -bPk-\W k nCi\) 



jeWi 



|VVf n w k \ + J2 7 iew* • - ' Wfc n 
jew, 



(29) 



Now, ignoring the content of the different sets, we can derive 
an expression for the expected edge set size of the primary 
graph. Note that the cardinality of the intersection of any 
two sets of given sizes, whose elements are unknown but are 
drawn from the same pool of elements, is a hypergeometric 
distributed random variable. Thus, we get: 

M 



E[|£„|] = J^EpA,] 



i=l 



^ M M ( 

2EE^ E [i w ' nWfe i] + ^ E 



i=i fc=i i 



E 

jeWi 



-E 



E J J-ew* ' ' Wfe n C * 
jew, 



1 M A/ f 



1=1 k=l 

k^i 



I V^i^fc i>k^iQk 
1 X X 



£ E[j j6Wjk .|w fc nA|] 

jew, 

(30) 



Note that the indicator function in the last term expectation 
can be only zero or one. Consequently, the expectation of 
its multiplication with |V\4 fl Ci\ can be only evaluated for 
Ij£TZ k = 1- If this is the case, this means that packet j 
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cannot be in the intersection of Wfe and £j. Consequently, 
this intersection is possible only with the other <pi — 1 packets 
and from the set of the remaining N — 1 packets. Thus, we 
get: 

N-l 

e [i jeRk ■ \w k n A|] = E n p [ J i6»* = x ' l Wfc n £ <l = "1 

2V-1 

= ^ n P \\W k n A| = n\l je n h = l] • P [/ ieW » = 1] 



71=1 



E 



n- 



Substituting (31) in (30), we get: 

M M 



(31) 



2 EE 

i=i k=i 



ipiipk , i>k^%Qk ^iQki>ki}Pi ~ 1) 



N 



N(N - 1) 









\ - V>/s 

TV 




fc=i 


A/ 


TV 




fe=i 



i + i - 



yi - l 
iV-1 



QkQi 
N-l 



(32) 



The first term in (34) is obviously the expected vertex degree 
of receiver i at time t. Now, we can derive the expected degree 
of receiver i after serving the maximal clique n as follows: 



E 



A (t+1) 1 — Ey / 
E 



E 



A (t+i) 



X 



fe=1 AT(iV-l) 



ip J^fc , (jffc(g fc-^fc) + -yg)(gi + -yo' 



+ E E * 



fcer„(ft) 

k^i 



E 



AT ' N(N-1) 
ipkX k (g. t + Xi) 



N(N — 1) 



a/ 



*;=i 



feGTp(K) 



fcer„(K) 



The expression for E 
approach. 



At+i) 



(36) 

can be derived using the same 



Appendix F 
Proof of Theorem 6 



Appendix E 
Proof of Theorem 5 

Using (25) in Theorem 6, we can derive an expression for 
the expected primary degree of receiver i G T(k) at time t+1, 
conditioned on the random vector X defined in Appendix A, 
as follows: 



E 



r A (*+i 



7» 



V>fc - X k ^ | (g* + X fe ) (ft + Xi) 



E 



fcer^(K) 



E %(i 



AT 

fc£T(«) 

Re-arranging the terms, we get: 



A' - 1 

(gfc + Xfc)(g i +X i ) 
AT - 1 

g fc (gi + jQ) 
AT-1 



(33) 



E 



A 



(t+i) 



A/ 



x = 



fc=l 



^ AT I N-l) 4-> 



QkQi 



M 



ipkQkX l 
^N(N-l) 

k^i 

(34) 



E £(i 



N - 1 



E 

k^i 



iPkXkjg. + X,) 
N(N - 1) 



(35) 



Consider an arbitrary vertex vij in the graph. From the 
adjacency conditions CI and C2 in Section IV, we can 
conclude the following facts: 

• Vertex v^ is not connected to any vertex of the same 
receive i. 

• If j £ Wfe, cannot be adjacent to any primary vertex 
of receiver k due to violation of C2, except for vertex 
Vkj which arises from CI. 

• If j € IZk, can be connected to any primary vertex 
of receiver k (induced from Wfe), except for all vertices 
Vki for which I ^ TZi I € 

From these facts, we can express the primary degree of a 



vertex Vij as follows: 



M 

E 

k=l 

kyii 



/ je w fc +W(|Wfe|-|>VfenA|) 



(37) 



where I x is an indicator function, which is equal to one if x 
is true and zero otherwise. 

Now, ignoring the content of the different sets, we can derive 
the expression for the expected primary degree of a vertex of 
receiver i. Consequently, we get: 



E [Ai] = E [Ay] 
E [IjeWk 



M r 

E 

k=l 
k^i 

M 



} + e [i jenk ] \w k \ - e [i je1tk ■ \w k n AO 



= E 

k=l 
k^i 



ipk , Qk4>k 



N 



N 



-E[l je1ik ■ \W k nd 



(38) 
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Substituting (31) in (38) and re-arranging, we get: 



A I 



k=l 
k=£i 

M 

E 

k=l 



tpk 
N 



1 + Qk 1 



QkQi 



Pi - 1 
N-l 



N - 1 



(39) 
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